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ABSTRACT 

Context. The search for radio emission from extra-solar planets has so far been unsuccessful. Much of the effort in 
modelling the predicted emission has been based on the analogy with the well-known emission from Jupiter. Unlike 
Jupiter, however, many of the targets of these radio searches are so close to their parent stars that they may well lie 
inside the stellar magnetosphere. 

Aims. For these close-in planets we determine which physical processes dominate the radio emission and compare our 
results to those for large-orbit planets that are immersed in the stellar wind. 

Methods. We have modelled the reconnection of the stellar and planetary magnetic fields. We calculate the extent of 
the planetary magnetosphere if it is in pressure balance with its surroundings and determine the conditions under 
which reconnection of the stellar and planetary magnetic fields could provide the accelerated electrons necessary for 
the predicted radio emission. 

Results. We show that received radio fluxes of tens of mjy are possible for exoplanets in the solar neighbourhood that 
are close to their parent stars if their stars have surface field strengths above 1-10G. We show that for these close-in 
planets, the power of the radio emission depends principally on the ratio (N c / B^ 3 ) 2 where N c is the density at the 
base of the stellar corona, and is the stellar surface magnetic field strength. 

Conclusions. Radio emission is most likely to be detected from planets around stars with high-density coronae, which 
are therefore likely to be bright X-ray sources. The dependence of stellar coronal density on stellar rotation rate and 
effective temperature is crucial in predicting radio fluxes from exoplanets. 

Key words, exoplanets, planetary radio emission, hot Jupiters 
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O ■ Searches for the radio signatures of extra-solar planets review f the c urren t status of theory and observation is 

j* ■ have not yet yielded any detections. The motivation for given in lZarkal (|2007f) 
. ' these searches was the well-known decametre emission from 

K> Jupiter which can outshine the radio emission from the One of the most promising reasons for this lack of de- 

$_| ' quiet S un at these wa velengths by some four orders of mag- tections is that tides raised on planets orbiting close to 

C3 ; nitude (jFarrell et all I1999D . This emission is believed to their host stars rapidly synchronise the axial spin with 

originate from non-thermal electrons travelling do wn the the orbital motion. In this case, the planet may be ro- 

converging f i eld lin es at the planet's magnetic poles (|Dulk|, tating slowly, producing only a weak magnetic field. Since 

|1985|; |Zarka|,|l998D. There are several possible sources for the emission frequency is proportional to the magnetic 

these electrons, including currents generated by corotation fi e id strength, this could result in emission at unobserv- 

breakdown inside the middle magnetosphere, the Io-Jupiter a hl y l ow frequencies. The other synchronisation that can 

interaction and the reconnection of Jupiter's magnet i c field take place is between the spin of the star and the orbit 

with the so l ar win d magnetic field (|Cowley fc Buncej|2001|; f the planet. This timescale is set by torques due to the 

|Saur et al.|, |2004[). Here we focus on the third possibility tidal b ulges raised on the star by the planet (see, e.g. lZahnl 

since it is common to all the magnetised planets. If these (| 1 977T ) and is generally expected to be longer than the 

electrons are accelerated to energies greater than 1-10 keV, sta r's main-sequence lifetime even for close-orbiting gas- 

they can be unstable to the electron cyclotron maser in- g i ant planets. The observed spin rates of planet-host stars 

stability and as a result they emit at the local electron are consistent with this expectation. A notable exception 

gyrofrequency /[MHz] = 2.8B(G). The observed peak of i s the late F star tau Boo, which rotates synchronously 

Jupiter's emission at 39.5MHz suggests a field strength of w } tri the 3.3-day orbit of its planet, whose mass is suffi- 

14.5G at the site o f the emission close to Jupiter's surface c i en t to lock the star ' s rotation within its mai n-sequence 

flConnerney et al.|, |1998|). The proposal that such a situ- lifetim e dMarcv et all. Il997t iLubow etall [l997h . More re- 

ation might occur in some of the known "Hot Jupiters" centlv iDobbs-Dixon et all (|2004f ) have considered the com- 

bined effects of tidal torques and angular momentum loss 

Send offprint requests to: M. Jardine from the star in the form of a wind. They suggest that syn- 



1. Introduction led to several searc hes but none of thes e has revealed the 

expected emission ([Rvabov et all [2004) . A comprehensive 
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Fig. 1. When the planet is far from the star, the stellar 
wind elongates the planetary magnetosphere in the direc- 
tion away from the star. R m is the extent of the plan- 
etary magnetosphere which is determine by the balance 
of pressures between the planetary and stellar plasmas. 
The shapes of the stellar and planetary magnetosphere are 
shown as dashed lines (not to scale). 



chronisation is more likely among F stars than among the 
lower mass G and K stars. 

The nature of the stellar wind may also have an im- 
pact on the expected radio emission, since the wind may 
compress the stellar magnetosphere a nd also may strip 
particles out of the planetary exosphere ( Gricfi meier et al. , 
2004 120051: 1 Jaritz et all 12005^ lLip~atov et al.l . l2005tlStevens1 . 
20051) . Fluctuations in the stellar wind pressure and speed 
that occur as a result of the stellar equivalent of solar coro- 
nal mass ejecti ons may, however, enhan ce any planetary 
radio emission (|Khodachenko et all l2007t ) . A recent search 
for radio emission from systems thought to be embedded 
in strong stel lar wind s has yielded no detections, but tight 
upper limits (jStevensl . 12003 iGeorge fc Stevens! . I2007D . 

It is clear from this modeling that conditions both exte- 
rior to the planet (such as the speed, density and magnetic 
energy of the stellar wind) and within the planet itself (such 
as the nature of the dynamo-generated magnetic field) may 
have an impact on the detectability of exoplanetary radio 
emission. There is another factor, however, which is the 
process by which the power carried in the stellar wind is 
converted into radio emission. Within our own solar sys- 
tem, there is a clear scaling between the magnetic or kinetic 
power in the solar wind i ncident on a planet and the ob- 
served output radi o power (|Farrell et al.l . ll999tlZarka et all 
120011; IZarkal 120071 ). This "radio-magnetic Bode's law" can 
be used to predict the radio emission from a variety of stel- 
lar wind or planetary magnetosphere conditions. 

In this paper we take a complementary approach. 
Rather than focussing on the physics of the stellar wind 
or the planet, we choose simple models for these and focus 
instead on the process of energy conversion. Our aim is to 
examine how the radio power scales with different stellar 
parameters, and to determine if this scaling is the same for 
all exoplanets. 

2. The model 

We begin by considering the interaction between the planet 
and its immediate environment. In the case of the outer 
planets, this takes the form of a stellar wind impacting the 
planetary magnetosphere. This wind carries both magnetic 
and kinetic energy, but it is the magnetic component that 
we are particularly int e rested in. Fo llowing previous work 
(e.g. lZarka eFaT]((200lT ); [Zarkal(l2Ml '). we can calculate the 
rate at which the stellar wind delivers magnetic energy to 
the planetary magnetosphere by calculating the Poynting 
flux (E x B//i) that is incident on the cross-sectional area A 




Fig. 2. When the planet is within the stellar magneto- 
sphere, the effective headwind due to the relative motion 
between the planet and the stellar magnetosphere elongates 
the plan etary magnetosph ere in the azimuthal direction 
(see also lZarka et al.l (|200lD ). 



of the planetary magnetosphere. We can approximate this 
maximum available power as P = AvB / fi. This depends 
on the cross-section of the planetary magnetosphere A = 
TtRm the wind velocity v and the stellar field strength B as 
measured at the planetary magnetosphere. This maximum 
available power sets an upper limit on the power that could 
be detected in radio emission. 

The situation for many "Hot Jupiters" is, however, 
somewhat different. For these planets, the orbital radius 
is extremely small - so much so that the planet may in 
fact be inside the stellar magnetosphere, where the stel- 
lar magnetic field lines are closed and confine the coronal 
gas. The maximum extent of the stellar magnetosphere is 
expected to increase with stellar activity (and hence with 
stellar rotation rate). Thus fr om a value of 2.5R(T) seen 
in eclipse images of the Sun jAltschuler &: Newkirk. Jr. I . 
119691: iBadalian fe Livshiti .ll986 : ). it may increase to ~19R0 
or 0.09AU for a star with rotation period about 6 days 
that is 10 times more active than the Sun in terms of 
surface magnetic flux density or soft X-ray luminosity 
(|Schriiver fc Title! . 120011: iSchriiver et all I2003T) . Given the 
youth of some exoplanetary systems, many planet host stars 
may have magnetic activity strong enough to support mag- 
netospheres extending to 0.1 AU or so. Hot Jupiters with 
orbital distance of order 0.05 AU will thus orbit within 
their host stars' magnetospheres. In this case, the velocity 
with which new magnetic field lines are carried into the in- 
teraction region is the relative velocity between the stellar 
and planetary magnetospheres. In the frame of the planet, 
however, this would appear to be similar to the effect of a 
stellar wind. 

The stellar wind therefore has two main effects on the 
planetary magnetosphere. The wind ram pressure com- 
presses the planetary magnetosphere and elongates it as 
shown in Figs |T]) and . The magnetic field that is carried 
in the wind is also pushed against the planetary magnetic 
field forming a narrow interaction region or current sheet. 
The magnetic field lines will pile up in this region until 
reconnection allows the field to dissipate. This process is 
observed both on the dayside and nightside of the Earth's 
magnetosphere and is related to the occurrence of aurorae 
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Fig. 3. Density of the gas in the region exterior to the 
magnetospheres of known exoplanets using orbital radii 
and stellar masse s taken from Th e Extrasolar Planets 
Encyclopaedia at htt p://exoplanet.eu/index.php| and as- 
suming R+ on M+. Close to the star, the planet is immersed 
in the the stellar magnetosphere, while far from the star, 
it is immersed in the stellar wind. At intermediate radii, 
planets may be in either regime and so no cases are plot- 
ted. 



Fig. 4. Magnetospheric radius (in Jupiter radii) for known 
exoplanets as a function of orbital separation. A stellar field 
B+ = 1G is assumed. 



3. Magnetospheric radius 

The first step in calculating the radio power emitted by 
exoplanets is to calculate the cross-section of the planetary 
magnetosphere that intercepts the stellar wind. The magne- 
tospheric radius is determined by pressure balance between 
the external medium and the planetary magnetosphere: 



Bp = Bl+np e v 2 . 



(1) 



as accelerated electrons stream along the Earth's magnetic 
field lines. 

The interaction of the planetary and stellar magnetic 
fields also generates an electric field at the interaction (or 
reconnection) site. If this electric field has a component 
along the magnetic field lines, it will accelerate electrons 
along the magnetic field lines. These electrons will mostly 
be slowed down by collisions, but some fraction may escape 
and be freely accelerated. It is this small fraction of the 
electron population that may provide the pool of energetic 
electrons needed to supply the electron cyclotron maser in- 
stability that would result in observable radio emission. 

In those rare cases, such as the anomalously massive 
t Boo b, where the stellar rotation is synchronised with 
the orbital motion, then there will be no relative motion 
between the stellar field and the planetary field and so al- 
though reconnection may happen sporadically due to the 
small velocity differences due to stellar differential rota- 
tion o r resistive instabilities such as perha ps the tearing 
mode (jFurth et all Il963t ICuntz et all . l2000h . it is unlikely 
to happen in the quasi-continuous manner that would be 
needed to power an observable radio signal. In the major- 
ity of observed cases, however, the star's rotation is not 
synchronous. New magnetic flux will then be brought into 
the interaction region at a rate determined by the velocity 
difference between the stellar and planetary fields. In this 
case the geometry of the interaction is altered and the elon- 
gation of the planetary magnetosphere will not be directed 
away from the star, but will trail the orbital motion ( see 
Figs. [2 and H]). 



here B p is the magnetic field of the planet which we take 
to be 14. 5G. We assume that this is a dipole, and so for a 
planet of radius R p this scales with distance as (R/R p )~ 3 . 
The field in the external region B e is the stellar field. Close 
to the star, where the field lines are still closed, we assume 
that this field is also a dipole and so for a star of radius 
R+ this field scales with distance as (R/R+)~ 3 . At distance 
greater than that of the stellar magnetosphere, however, 
where the field lines have been opened up by the stellar 
wind, the field scales as (-R/i?*) -2 . For the inner planets, 
the velocity in ([1]) is the relative velocity between the star 
and planet: 



3 x 10 4 



0.03AU 



lid 



P* 

Porb 



where we take the stellar r otation period to b e 11 days 
(appropriate for HD189733 (|Bouchv et all I2005D ). For the 
outer planets, this velocity is the wind velocity. This is an 
extremely difficult parameter to determine observationally, 
though use of hydrogen-wall absorption to de duce the wind 
ram pressure promises to be a useful tool (|Wood et all 
120041 ). Models of the s e winds also give very confl icting re- 
sults (|Stevend 120051 : iHolzwarth fc Jardinel 12007ft and so, 
for the moment, we simply assume that the winds have 
reached their terminal velocity and so the wind density falls 
off with distance as (i?/i?*) -2 . We scale the stellar wind ve- 
locity to solar values, assuming a nominal solar wind speed 
of SOOkms- 1 and a density at 1AU of 1.7x 10~ 20 kgm- 3 
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(|Cox et al.Lfl99lh . This gives a number density 



N = N W (1AU) 



= 10' 
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1AU 



1 -2 



(2) 
(3) 



For close-orbiting planets, however, the external density 
is not that of a wind, but that of the stellar magnetosphere. 
This is likely to be a function not only of the orbital radius 
but also the magnetospheric temperature and the stellar ro- 
tation rate (which influences the level of magnetic activity) . 
We assume that the stellar magnetosphere is in isothermal, 
hydrostatic balance and so obtain (see Fig. [3]) 



N = TVoexp 



/iihh GM* 
kT R* 



1 



R* 
R 



orb 



(4) 



with the number density at the base of the corona N = 
4 x 10 14 m~ 3 , the tem perature T=1.4 x 10 6 K and the mean 
particle mass [i = 0.6 (jBadalian fc Livsh its. 1986). We note 
that this gives a different scaling of the density N with 
orbital radius in the two cases; for the outer planets, N oc 
Br 2 , whereas for the inner planets, N cx R . 

Usin g these figures, we find, in com mon with pr evious 
studies ijCuntz et all 120001: IZarka et all . 120011: IZarkal I2007D 
that for the more distant planets, the contribution of the 
external field to pressure balance is negligible. For the inner 
planets the opposite is true, and we may neglect the ram 
pressure associated with the relative motion of the planet 
and the star. This then gives, for the inner planets: 



R„ 
R,, 



15.7 



Ro 



0.03AU 



B n 



-i i/3 r 



14. 5G 



R, 



Rr. 



-B* 
1G 



-1/3 



and for the outer planets: 



Rm 

Rp 



46 



Rorb 


1/3 


Bp 


1/3 




5AU_ 




14.5G_ 




-500fcms _1 . 



-1/3 



(5) 



(6) 



We show in Fig. (T4|) the corresponding sizes of the planetary 
magnetospheres (assuming that all planets have a magnetic 
field similar to Jupiter's). While for planets far from their 
parent stars the magnetospheric size is independent of the 
the stellar field strength (being determined mainly by the 
wind ram pressure), for the close-in planets, the pressure 
of the stellar magnetic field can have a significant effect 
in compressing the planetary magnetosphere. If the stel- 
lar field strength is too high, the planetary magnetosphere 
may be totally crushed, but for intermediate stellar field 
strengths or stronger planetary fields the planet may still 
retain a magnetosphere even in a very close orbit. 



4. The magnetic interaction between the star and 
the planet 

The interaction of the planetary magnetosphere with the 
external medium not only determines the size of the plan- 
etary magnetosphere but it also generates an electric field 
at the site where the magnetic fields of the planet and the 
star meet given by 



E = -v x B +j/H 



(7) 
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Fig. 5. The functional form of the dependence of the num- 
ber of runaway electrons on the ratio of the applied electric 
field to the Dreicer field. 
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Fig. 6. Ratio of the electric field E developed at the re- 
connection site to the Dreicer field Ed. Electron runaway 
occurs when E > Ed- A stellar field B+ — 1G is assumed. 



where the current density j = V x B //i. Here we have scaled 
the electric field to its value far from the interaction site 
Eq = —voBo where vq and Bq are the velocity and magnetic 
field strength. The current density is scaled to Bo / I^L where 
L is the typical lengthscale for the interaction which we 
take to be 0.1i?j. The results are in fact fairly insensitive 
to the exact value of L chosen. The magnetic Reynolds 
number H = vqL/tj where the Spitzer diffusivity is given 
by r) = IOST^/Ws- 1 (jPriestl . [l98l . 

Thus, well away from the interaction site, the electric 
field is simply determined by the velocity and magnetic 
field strength and we have E = — v x B, while at the centre 
of the interaction region where v — 0, the electric field is 
determined by the magnetic Reynolds number and we have 
E = j/#. 

It can be shown (see Appendix for details) that there 
is a component of this electric field that is parallel to the 
magnetic field and so is capable of accelerating electrons 
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Fig. 7. Power in electrons accelerated by the interaction 
of the stellar and planetary magnetic fields. A stellar field 
-B* = 1G is assumed. 
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Fig. 8. The minimum stellar magnetic field needed to accel- 
erate electrons to above lOkeV in the interaction between 
the magnetic fields of the star and planet. 



along the magnetic field lines. As shown in the Appendix, 
this electric field (En) has a maximum at the interaction 
site between the fields of the star and the planet. Its mag- 
nitude falls away with increasing distance z from this site 



as 



E 



(8) 



H = E B = E\\{Q)e~ z H . 

It is the magnitude of this parallel electric field that is im- 
portant in accelerating electrons in the interaction region 
Its value is given simply by 



£n(0) = - 



2V2vr 



(9) 



r(i/4) 

Hence, although there is a weak dependence on the mag- 
netic Reynolds number (and hence on the temperature of 
the local plasma) , the magnitude of this accelerating elec- 
tric field is determined primarily by the magnetic field 
strength and the velocity at some distance from the in- 
teraction site. 

5. Power in accelerated electrons 

Once the electrons are accelerated by this field however, 
they may be slowed down by collisions with a collision fre- 
quency 

v c = 7.5 x l(T 5 T- 3/2 n (10) 

that depends on the local temperature T and d ensity 
n, here given in SI units dKruskal fc Bernstein! . 1 1964 



Tandbcrg -Hanssen fc Emslid . ll988T ). However, if the electric 
field is stronger than a critical value known as the Dreicer 
field Er>, electrons may be freely accelerated and so will 
"run away" out of the thermal distribution. The Dreicer 
field is also a function of the local temperature and density 
and is given by 

E D = 18 x lQ- 12 nT- l \ymT x ]. (11) 
The number density of runaway electrons is then given by 
iV run = 0.35nvJ{E/E D )[m- 3 ] (12) 

= 2.6 x io-Vr- 3/2 /(£/£D)[ur 3 ] (13) 



where the function /(E/Ed), shown in Fig. [5J is given by 



(_E_\ = (Ed 
\E D ) I E 



3/8 



exp 



( 2E D 
V E 



1/2 



Ed 
~\E 



(14) 



The behaviour of this function governs the way in which 
the electric field generated at the interaction site can affect 
the number density of electrons that run away and hence 
can power the radio emission. Once this imposed electric 
field is greater than the Dreicer field, there is a steep rise 
in J(E/Ed) with increasing E. Once the ratio E/Ed is 
sufficiently large, however, / becomes insensitive to further 
increases in E/Ed and saturates at a value of about 0.3. 
We can see from Eqs . l9l and ITTI t hat for a given temperature, 



E v B 
oc . 

E D n 



(15) 



The local temperature at the acceleration site n will depend 
on the details of the reconnection process and its determi- 
nation would require a full analysis of the local dynamics. 
It is a reasonable assumption however that it scales with 
the density of the plasma flowing into the reconnection site 
such that n oc N. We determine the constant of propor- 
tionality empirically by fitting the predicted radio power to 
that observed for the solar system planets, and find it to 
be of order 7. Hence, as shown in Fig [HI the ratio E/Ed 
is the same for all the solar system planets since the solar 
wind velocity is a constant and both the density and field 
strength vary as R~^,- This is not the case for close-in exo- 
planets, where the velocity, density and field strength vary 
differently with orbital radius. 

Electrons that are accelerated over a lengthscale R m to 
an energy K = qER m carry a power 



P = nR 2 m vN run K[W} 



(16) 



where the velocity v is either the wind velocity (for the outer 
planets) or the relative velocity of the star and planet (for 
the inner planets). Hence P oc w7V 2 i? 2 n and so the domi- 
nant factors determining the output power are the size of 
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the planetary magnetosphere and the density. Using expres- 
sions [3] to |51 we can show that for given planetary parame- 
ters, this power scales as 



P cx N 2 vB± 2/3 



for the inner planets and 



P cx NlvRtR-JI* 



(17) 



(18) 



for the outer planets. Both of these cases are shown in 
Fig. \7\ For those planets embedded in the stellar wind, 
the fall-off in power with increasing orbital radius is very 
pronounced. For the inner planets, however, the power is 
almost independent of orbital radius. As the orbital radius 
decreases, the number of electrons available to be accel- 
erated rises (since N run cx n 2 cx R~ r \)- This is balanced 
however by the shrinkage of the collecting region of the 
magnetosphere that is apparent in Fig. [4] since cx R 2 rh - 

The power carried by the accelerated electrons for these 
inner planets is therefore largely independent of the plan- 
etary orbital radius (albeit with a weak variation due to 
the relative velocity v). The magnitude of this power is 
however affected by the stellar parameters we have chosen. 
For a stellar field strength of only 1G, not only is the typ- 
ical electron energy of the runaway electrons reduced to 
a few keV, (and so we set the characteristic electron en- 
ergy to leV) but also since the electric field generated is in 
many cases less than the Dreicer field, fewer electrons are 
accelerated. A higher stellar field strength would of course 
provide a larger electric field and hence accelerate electrons 
to higher energies, but it would lead to a greater compres- 
sion of the planetary magnetosphere and hence reduce the 
upper limit for the lengthscale over which these electrons 
could be accelerated. 

We can quantify the effect of the magnetic field by deter- 
mining the minimum stellar field strength needed to ensure 
that the distance L min required for electrons to be acceler- 
ated above (say) lOkeV is less than the maximum available 
lengthscale (R m ). Thus 



R 



oc 



3/2 
orb 



' mm ~ R 3 

for close-in planets, while for the outer planets 



R 



5/3 



(X 



m 



(19) 



(20) 



We show this minimum stellar field strength in Fig. ([8|) 
where we assume a stellar orbital period equal to that of 
HD 189733 for close-in planets and a stellar wind speed of 
500kms _1 for the outer planets. Clearly, only a small in- 
crease in the stellar field strength above 1G would be ade- 
quate to ensure that a pool of electrons of energies greater 
than lOkeV would be available to power the electron- 
cyclotron maser instability for most exoplanetary systems. 
Even below this limit, energies of greater than IkeV are 
easily achievable. 

6. Comparison with the solar system 

While Fig. [7J illustrates the scaling of the expected power 
with orbital radius, the level of this power depends on our 
assumptions for the stellar field strength and densities. The 
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Fig. 9. Radio power assuming the electron cyclotron maser 
instability converts power in accelerated electrons into radio 
power at 10% efficiency. Values for the sol ar system pla nets 
(shown as open squares) are taken from IZarkal (|2007f ). A 
best-fit line through these values is shown. A stellar field 
= 1G is assumed. 



fraction of this power that may be converted into observ- 
able radio emission also depends on the details of the elec- 
tron distribution and the elec t ron c yclotron maser insta- 
bility (see e.g. revi ew by (fWul. 1 19851) 1. Values of 1 - 20% 
have been claimed (lLouarn . Il992t iMacKinnon et all Il992t 



ICairns fc Binghaml 120021 IVorgul et all . I2005D and so we 
assume that 10% of the power in accelerated electrons is 
converted into radio power. This allows us to compare our 
predicted radio power with that observed from the planets 
in our own solar system, which is directly proportional to 
the input magnetic power carried by the solar wind (see 
Fig. ^j). This "radio-magnetic Bode's law" has a constant 
of proportionality or efficie ncy of 1 — 10 x 10 -3 relating the 
magnetic and radio powers (|Zarkal . 120071 ). We can write this 
efficiency e = P r adio/-Pma g netic as 



0.1 



B 2 



(21) 



where K is the typical energy acquired by the accelerated 
electrons. The factor of 10% is an estimate of the fraction 
of the power in accelerated electrons that can be converted 
to radio emission. Since we expect that N lun cx N 2 , this 
gives 



(22) 



For the outer planets (such as those in our own solar sys- 
tem) we can easily reproduce this constant efficiency e by 
noting that both TV and B scale as (R or b/Ri)~ 2 This gives a 
constant efficiency whose magnitude depends principally on 
the local field strength and density. If we scale these to their 
values at the Earth's orbit (a density of 1.7 x 10~ 20 kgm -3 
and a field strength of 3.5nT) then we can reproduce the 
emission from the solar system planets shown in Fig. [9] 
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Fig. 10. Planetary flux at 40 MHz (in mJy) received at 
Earth. A stellar field = 1G is assumed. In order of 
increasing orbital radius, the systems that lie above lOmJy 
are: Gliese 876 d, HD 189733 b, 55 Cnc e, HD 179949 b, 
Tau Boo b, 51 Peg b, Ups And b, HD 69830 b, and HD 
160691 d. 



7. Discussion 

For the exoplanets that are far from their parent stars, a 
universal "radio-magnetic Bode's law" might therefore be 
expected, but with a slope that is determined by the ra- 
tio of N/ B in the stellar wind. These two factors may well 
not be independent and indeed may vary with both spec- 
tral type and stellar age. IGriefimeier et al.1 (|2005f ) have con- 
sidered the effect on the detectability of radio emission of 
changes in stellar winds as stars age. They conclude that 
young systems with more powerful winds are the most likely 
candidates. 

The inner planets are, however, much more promising 
candidates for detection because of their larger predicted 
radio powers. The magnitude of this power also depends 
on their coronal density and to a lesser extent the stellar 
field strength. Increasing B* increases the input power and 
so moves all the points in Fig. [5] to the right. The output 
power, however, is not significantly affected. Thus, for ex- 
ample, if we consider HD189733 we can conv ert the radio 
power to a flux j> received at Earth follow ing iFarrell et al.l 
(| 19991 ) and IGriefimeier et all (|2005l . l2007al ) to obtain 



4> 



nd 2 5f 



(23) 



Here d is the distance to the star and we assume planetary 
properties similar to Jupiter (fZarka, 2004]). Hence Sf is the 
bandwidth, taken as Sf — f c where f c — 40MHz is the max- 
imum cyclotron frequency and we assume that the emission 
is beamed into a solid angle O = 1.6 sr. While for a stellar 
field of 1G (as shown in Fig. [TDJ( , the flux from HD189733 
is only 14mJy, if we use the observed stellar field strength 
of 40G, the electric field increases so that the typical elec- 
tron energy rises, but the magnetosphere shrinks, with the 
result that the flux changes to only 15mJy. This assumes, 
however, that the coronal density is equal to that of the 
Sun. If we assume s calings in the range TV o c O 0-6 to N oc Q 
(|Unruh fc Jardind . Il997t llvanova fc Taaml I2003D then for 



the rotation rate of 11.73 days appropriate for HD189733, 
the flux rises to from 15mJy to 39 - 72 mJy respectively. 

This example points to the importance of the stellar 
coronal density in determining the possible radio power. It 
suggest that those stars with high X-ray flux (which like the 
radio flux also scale as N 2 ) are the most likely candidates. 
Although rapid rotation may therefore lead to a greater ra- 
dio power, most radial- velocity searches for exoplanets have 
tended to exclude rapidly-rotating objects from their survey 
samples, because their high rotational broadening degrades 
the precision with with radial velocities can be determined, 
and their surface activity causes radial- velocity "jitter". 
Stars that are magnetically active tend to have significant 
radio emission of their own, which c ould mask any plane- 
tary signature (jBenz fc Guedell . ll994D . Indeed, we might ex- 
pect planet-induced emission due to the electron-cyclotron 
maser process at the star. This would appear at the cy- 
clotron fr equency of t he stellar magnetic field f[MHz] — 
2.85* [G] ((Dulkl . H985T) . Th us, for HD189733 wi th an ob- 
served field strength of 40G (jMoutou et al.l . l2007t) we might 
expect emission at 112MHz. In order for this emission to 
escape from the stellar magnetosphere, however, the local 
cyclotro n frequency must be greater than the plasma fre- 
quency ( Dulk .ll985tlBingham fc Cairnsl l2000t IZarkal l2007t 
IGriefimeier et all . l2007bfh This means that at the site of 
emission, 



B>3x 10- 6 A e [m- 3 ] 1/2 [G]. 



(24) 



For emission from the planet, this condition is easily met, 
requiring only a planetary field strength of 0.1 G for a den- 
sity of 10 9 m -3 . If we are considering emission from the 
star, however, this condition requires that > 60G for a 
density of 4 x 10 14 m -3 . 

Stars that have higher field strengths will have a higher 
emission frequency and so may be easier to observe. In prin- 
ciple, this coul d be distinguished f rom the intrinsic emis- 
sion of the star (|Kellett et all [2002) by its modulation due 
to the orbit of the planet, but this would require obser- 
vation over many orbits to average out the effect of tran- 
sient stellar emission due to stellar coronal activity. This 
process may also lead to chromospheric signatures of the 
interaction as electrons from the reconnection site are ac- 
celerated a long the stellar field lines and im pact the stel- 
lar surface. IShkolnik et all (|2003l I2005L f2008f ) have reported 
plan et-induced chromosp heric activity on HD 179949 and v 
And. lMcIvor et all (|2006l) modelled this process and showed 
that this interaction could explain the chromospheric activ- 
ity, but not the phase lag betw een the planet a nd the chro - 
mospheric hotspot for v And. iPreusse et al.l (|2005l . [2006) 
have, however, reproduced this phase lag by considering 
the propagation of Alfven waves generated by the passage 
of an exoplanet through a stellar magnetic field. 

The level of magnetic activity may be important in de- 
termining the rate at which the stellar field changes may 
also be important. These changes are most likely to be 
driven by new flux emergence through the stellar surface, 
or by surface transport processes: diffusion, differential ro- 
tation and meridional flows. While there is now a general 
acceptance that surface differential r otation decreases with 
decreasing mass (|Barnes et all 12005b . the nature of the sur- 
face diffusion or meridional flows is as yet uncertain. 

The stellar magnetic field also has a direct effect on 
the emission process through its control of the size of 
the planetary magnetosphere. Indeed, if the stellar field 
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strength is too large, the planetary magnetosphere may be 
totally crushed. We can calculate the minimum stellar field 
strength required to do this by setting the size of the plan- 
etary magnetosphere equal to the planetary radius in the 
condition for pressure balance ([T]) , and setting the ram pres- 
sure term to zero. This gives 



B*,min = 3.8 x 10 3 



Rorb \ ' 

0.03AU J 



R* 



[G]. 



(25) 



For the subset of planets that transit their host stars, 
there is the possibility of observing radio occultations as 
the planet passes behind the star. With the observation 
of repeated occultations any time-dependence in the emis- 
sion due to variations in the emission process itself could 
be averaged out. Such observations might provide informa- 
tion about the structure and orientation of the planetary 
magnetic field, since most of the emission is expected to 
originate in the region of the magnetic pole close to the 
planetary surface. Even if this emission is present, how- 
ever, it may be beamed out of the line of sight. Recent 
claims of electron-cyclotron emission from low-mass stars 
with kG magnetic fields are, however, extremely encour- 
aging and suggest that occultations are indeed possible 
as the source region rotates out of view b ehind the star 
(Hal linan et "all 120081; lAntonova et aT1 . l2008h . The emission 
is at much higher frequencies than would be expected from 
a planetary magnetic field however. 



8. Conclusions 

In studying the possible radio emission from exoplanets, we 
have focussed on the electric field that is generated when 
the planetary and stellar magnetic fields reconnect. This 
electric field can accelerate electrons to energies sufficient to 
power the electron-cyclotron maser instability which is be- 
lieved to be responsible for much of the radio emission from 
the planets in our solar system. We have assumed that the 
mass, radius and magnetic field strength of all exoplanets 
is the same as Jupiter, but used the observed stellar values. 
By modeling the electric field developed in this interaction, 
we have calculated the number of accelerated electrons pro- 
duced. We find that it scales as the square of the density N. 
Thus the output radio power scales as N 2 while the input 
power from the stellar magnetic field scales as B 2 . Since 
for planets immersed in a stellar wind (as is the case for 
the planets in our own solar system) both the density and 
the field strength scale as R 2 rh , this naturally predicts that 
the output and input power should be proportional, as is 
observed for the solar system. 

We find, however, that the output power from the in- 
ner planets behaves differently. The planets lie inside the 
closed magnetic field (or magnetosphere) of the star. The 
magnetospheres of these planets are confined not by the 
ram pressure of the stellar wind, but by the magnetic pres- 
sure of the stellar magnetosphere. The density and mag- 
netic pressure within which they are immersed do not vary 
with orbital radius in the same way as for the outer planets, 
and so we would not expect the output and input powers 
to be proportional. In fact, we find that as the planetary 
orbital radius decreases, the density rises and so the pool 
of available electrons increases, but this is exactly balanced 
by the shrinking of the magnetosphere which reduces the 



cross-section of the planet. The result is that while the in- 
put power might rise, the output power saturates at a value 
determined by the ratio 



Pout oc 



1/3 



(26) 



where N c is the density at the base of the stellar corona and 
B+ is the surface field strength. We conclude, therefore, that 
the radio power emitted by "hot Jupiters" is determined 
by the way that coronal density scales with stellar field 
strength as a function of stellar rotation rate and effective 
temperature. 
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Appendix A: Modelling the reconnection of the 
stellar and planetary magnetic fields 

We model the interaction of the stellar and planetary mag- 
netic fields in t he manner similar to that used for colliding 
wind systems (pardine et all |l996). Because of the large 
length scales involved we neglect the radius of curvature 
of the reconnecting field lines. In this case, the equation of 
motion of the flow is just 



p(v • V)v = -V [p + B 2 j2p\ + pvV\, 



(A.l) 



where v is the kinematic viscosity and the other symbols 
have their usual meaning. For this local analysis, we use a 
Cartesian coordinate system where the z-axis lies along the 
line of centres of the star-planet system (see Figure IA. 1|) . 
As the magnetic fields of the star and planet are pushed 
together by the local flow, a current sheet is generated. 
Within this current sheet, components of the magnetic field 
that are antiparallel will annihilate, while components that 
are parallel will combine. We therefore chose to align the x- 
and y-axes such that the x-components of the magnetic field 
on either side of the current sheet are exactly antiparallel 
and hence cancel out within the current sheet, while the 
y-components are parallel. As a result, at the center of the 
current sheet the magnetic field has only a y-component. 
The magnetic field is then of the form 



consider only the case of a zero-vorticity flow for which an 
exact solution for the ma gnetic field profile can be found 
(jSonnerup fc Priesti Il975f h In this case, the velocity is of 
the form 



(x,y,-2z) 



(A.4) 



where all distances are scaled to some global lengthscale L 
so that x = xjh etc and v = v/vq where vo is the velocity 
with which field lines are being carried towards the current 
sheet. Now, by our assumption of a steady state, V x E = 
where the electric field is given by a simple Ohm's law 



E = v x B + j/cr. 

This can be written as 



B = 0, 



(A.5) 



(A.6) 



where H = vqL/t] is the magnetic Reynolds number and 
B = B/B a where B Q = B(z = 1). The electric field E = 
E/vqBq has components 



E = {-2zB y - B y /H, 2zB x + B x /H, -xB y + yB x ) 

If we take the component of E along B to be E 
straightforward to show that 



#||(0)e 



(A.7) 

, it is 

(A.8) 



E| | — E • B 

and so we would expect that the centre of the current sheet 
is the site of the most efficient particle acceleration. We 
therefore concentrate our attention on the particles accel- 
erated there. The degree to which the magnetic field rotates 
as it comes into the current sheet then depends on the rel- 
ative magnitudes of the x- and y-components outside the 
current sheet. We specify this as one of the boundary condi- 
tions, in fact making them equal at z — R m so that the field 
rotates through 7r/4 radians. With boundary conditions 

B x (0) = 0, B x (l) = l, By(0) = 0, 5 y (l) = l (A.9) 

the solution to (|A.6[) can be expressed in terms of modified 
Bessel functions as 



Bx = 



and 



z l/2 e (l-z 2 )H/2 

h/*(H/2) 



-I 1/4 (z 2 H/2) 



y l/2 e (l-z 2 )H/2 



-I_ 1/4 (z 2 H/2) 



For H » 1 we can write I±i/ 4 (H /2) ss (e 
that at the centre of the current sheet we have 



(A.10) 

(A.ll) 

H /nH) 1/2 so 



By(0) 



S2tt 



r(3/4) 



B = (B x (z),B v (z),0) (A.2) 

and the current density j = V x B/fi is 

M = (-B' y ,B' x ,0). (A.3) 

The variation with distance from the centre of the current 
sheet of the velocity and magnetic field for colliding mag- 
netofluids has been determined for arbitrary vorticities by 
iJardine et all (|1993f ).with a numerical solution that shows 
that the magnetic field strength at the centre of the cur- 
rent sheet varies only slowly with vorticity. We therefore 



B x (0) 

and 

E(0) = 
or 

25(0) = 



2V2tt 

rW4) 

2V2^ 

~rW4) 

2V2^ 



#3/4 



H 



-1/4 



H-^voBo. 



(A.12) 
(A.13) 

(A.14) 
(A.15) 
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Fig. A.l. Schematic diagram of the coordinate system used 
for the magnetic field at the centre of the current sheet 
where the stellar (grey) and planetary (black) magnetic 
fields reconnect. The z-axis points along the line of centers 
joining the star and planet. 



